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INTRODUCTION 


Satellite  to  sea  surface  altimeter  data  can  be  used  to  estimate  sur¬ 
face  gravity  (Ag)  in  ocean  areas  by  assuming  the  mean  sea  surface  (f) 
measured  with  the  altimeter,  suitably  corrected,  is  the  geoid  (N).  The 
Geos-3  and  Seasat-1  radar  altimeters  have  demonstrated  the  concept  and 
Geosat-A  will,  hopefully,  provide  systematic  and  complete  coverage  of  the 
world's  oceans.  We  will  discuss  the  procedures  and  errors  to  be  expected 
with  Geosat-A  in  three  broad  categories. 

1.  Measurement  of  Sea  Surface 

(a)  Accuracy  of  the  altimeter 

(b)  Accuracy  of  Geosat-A  ephemeris 

2.  Non-geodetic  contributions  to  observed  sea  surface 

(a)  Oceanographic  signals 

(b)  Environmental  corrections 

3.  Integral  for  calculations  of  Ag  from  N 

(a)  Integration  limited  to  a  spherical  cap 

(b)  Error  propagation 

In  general  these  are  all  complex  issues.  In  particular  2. a  and  2.b  involve 
much  of  what  is  known  as  physical  oceanography  and  tropical  meteorology. 

We  will  narrow  the  discussion  by  focusing  on  a  particular  objective  and 
data  analysis  strategy;  viz.  calculation  of  Ag  in  the  world's  oceans.  One 
can  describe  this  as  follows.  The  end  product  is  a  global  map  of  gravity 
at  sea  with  a  resolution  of,  say,  15  minutes  (15').  This  gravity  map  will 
be  determined  from  a  geoid  surface  obtained  from  combining  data,  taken 
over  more  than  10  years,  from  Geos-3,  Seasat-1,  and  Geosat-A.  The  compu- 
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tation  will  be  made  with  a  Molodensky  type  integral  using  data  within  a 
spherical  cap  of  a  few  degrees.  The  computation  will  use  a  reference  sur¬ 
face,  say  WGS83,  truncated  at  some  degree  and  order,  and  aim  for  an  accuracy 
of  a  mi 1 1 i gal  or  better. 

The  accuracy  of  recovered  gravity  anomalies  depends  on  the  accuracy  of 
the  geoid  and  the  resolution  sought  as  illustrated  in  Table  1.  Table  1  il¬ 
lustrates  the  relation  between  geoid  and  gravity  accuracy  for  a  number  of 
block  dimensions  a  and  b.  Since  the  relation  scales  it  can  be  used  for  any 
assumed  accuracy.  For  example,  for  a  10km  x  10km  square,  a  geoid  uncer¬ 
tainty  of  M=lcm  would  produce  a  gravity  anomaly  map  of  n=0.76mgal. 

The  mathematical  relation  to  calculate  a  gravity  anomaly  (Ag)  from 
geoid  height  (N)  is  a  (singular)  integral  equation  due  to  Molodensky  (Molo¬ 
densky,  et  al.,  1962).  The  error  in  this  estimate  has  two  sources:  First, 
observation  errors  (e)  and  real  physical  departures  of  the  geoid  (N)  from 
the  mean  sea  surface  (  f ■ ) ,  and  second,  errors  in  the  numerical  implementa¬ 
tion  of  the  integral.  The  first  error  source  is  due  to,  for  example, 
satellite  orbit  error,  ocean  currents,  tides,  storm  surges,  salinity  and 
temperature  variations,  and  eddies.  In  general  these  vary  with  time  and 
position,  and  have  a  maximum  amplitude  of  about  one  meter.  The  large  am¬ 
plitudes  have  time  variability  of  weeks  to  months,  and  with  data  taken  for 
several  years,  they  may  average  out.  These  'errors'  will  be  discussed  in 
this  Report. 

The  second  error  source  has  two  parts.  The  first  is  due  to  the  fact  that 
the  area  of  integration  for  the  Molodensky  integral  is  the  whole  unit 
sphere.  In  other  words  this  calculation  assumes  complete  coverage  of  the 
geoid.  Forming  the  complete  integral  is  impractical.  In  addition,  the 


altimeter  only  provides  geoid  data  in  the  oceans.  Therefore,  the  integra¬ 
tion  is  limited  to  a  spherical  cap  around  the  integration  point.  The  error 
due  to  neglecting  this  outer  zone  can  be  quantified.  Secondly,  the  integral 
must  in  practice  be  replaced  by  a  finite  summation.  This  error  of  quantiza¬ 
tion  must  also  be  considered.  However,  the  integration  error  from  limiting 
the  area  of  integration  to  a  spherical  cap  can  be  made  less  than  a  mi 11 i gal 
by  use  of  a  suitable  reference  geoid  and  spherical  cap  size  (Sections  II, 

IV,  V,  and  VI;  Gaposchkin,  1983).  We  proceed  to  analyze  the  Molodensky 
integral,  and  examine  an  implementation  of  this  integral. 
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ESTIMATE  OF  ACCURACY  N  — >  dg 


We  begin  by  finding  a  general  relation  between  the  accuracy  of  mean 
gravity  anomalies  (Ag)  and  mean  geoid  heights  (N)  for  areas  with  sides  a 
and  b  assuming  complete  knowledge  of  N  or  Ag.  Following  Moritz  (1974)  we 
start  with  Stokes'  formula: 


N  = 


4  Try 


t 


S  (*)  Ag  dor 


(1) 


where  R  is  the  earth's  mean  radius,  y  is  the  mean  value  of  gravity,  S  (^) 
is  the  Stokes  function,  and  a  is  the  unit  sphere.  We  set  by  error  propaga¬ 
tion: 


M2  = 


8*-  y 


-g  J [S(*)]2  sin  ^  d  ^ 


(2) 


where  M  is  the  standard  error  of  the  geoid  height  (N).  It  is  shown  (Moritz, 
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1974)  that  for  square  blocks  of  side  a,  with  L=a  m  , 


a  2  2 

g  m 
4?r  yt 


*  <"o> 


(3) 


where  M  is  the  standard  error  of  the  gravity  anomaly  (Ag),  and 


We  generalize  (3)  and  (5)  for  rectangular  block  -'f  sides  a  and  b  as: 


M2 


a  b  m 


4ir  7 


*(n0) 


(6) 


y  R 

nQ  =  \/ab~ 


(7) 


The  function  *(nQ)  is  very  slowly  changing  and  in  our  range  of  interest 
is  approximately  20.  Equation  (6)  can  be  used  to  estimate  the  gravity  ac 
curacy  given  a  geoid  accuracy,  or  to  estimate  a  geoid  accuracy  given  a 
gravity  accuracy.  Table  1  lists  values  for  various  M,  m,  a,  and  b. 


TABLE  1 


UNCERTAINTY 

FOR  GRAVITY 

ANOMALY  VS.  GEOID 

HEIGHT 

a 

b 

no 

M 

m 

1M 

1M 

1c ml 

(mqal ) 

1 

1 

20000 

25.173 

10 

69.94 

10 

10 

2003 

21.292 

13 

10.0 

10 

10 

2003 

21.292 

10 

7.6 

25 

12.5 

1133 

20.151 

10 

4.55 

20 

20 

1001 

19.903 

26 

10.0 

20 

20 

1001 

19.903 

10 

3.93 

25 

25 

801 

19.456 

31 

10.0 

25 

25 

801 

19.456 

10 

3.18 

50 

50 

401 

18.067 

61 

10.0 

50 

50 

401 

18.067 

10 

1.65 

100 

100 

200 

16.665 

118 

10.0 

100 

100 

200 

16.665 

10 

0.86 

For  example,  a  15'xl5‘  grid  at  a  latitude  of  60  degrees  would  be  a  rectangl 
of  approximately  28km  x  18km.  Ten  cm  geoid  heights  would  correspond  to 
4.55  mi 11 i gal  gravity  anomalies.  Since  these  variances  scale,  1  cm  geoid 
heights  would  correspond  to  0.455  mgal  gravity  anomalies,  and  25cm  geoid 
heights  would  correspond  to  11.4  mgal  gravity  anomalies.  These  estimates 
are  confirmed  very  well  in  the  simulations  using  the  Molodensky  integral. 


MODIFIED  MOLODENSKY  INTEGRAL  AND  ITS  ACCURACY 


The  Molodensky  integral  (Molodensky,  et  al . ,  1962)  is  generally 
written  as: 


where  Agp  and  Np  are  the  gravity  anomaly  and  geoid  height  at  the  computa¬ 
tion  point  p  and  rp  is  the  cord  distance  between  the  point  p  and  the  dif¬ 
ferential  area  element  d<r  .  Now: 

rp  =  2R  Sin  (  */2)  (9) 

where  ^  is  the  central  angle  between  point  p  and  the  area  element  dot  ,  and: 

2 

dor=  R  Sin^d^da 

The  Molodensky  integral  is  a  global  integral  and  in  principle  requires 
knowledge  of  N  everywhere.  In  addition,  in  practice  the  integral  (8)  is 
replaced  with  finite  sums.  Also,  as  a  practical  matter  N  is  not  known 
everywhere  and  considerable  simplification  and  economy  is  achieved  by 
limiting  the  summation  to  geographical  regions  near  the  computation  point 
p.  The  quantization  error  is  explored  here  with  a  numerical  simulation. 


The  effects  of  neglected  zones  outside  a  spherical  cap  has  been  in¬ 
vestigated  by  Buglia  (1976).  Buglia  proceeded  to  write  the  integral  (8) 
in  two  parts:  one  for  the  integral  within  a  spherical  cap  of  angular  radius 
and  a  second  part  for  the  remaining  area  of  the  unit  sphere,  which  we 
call  the  outer  zone  (Buglia,  1976,  eq  3).  He  calls  the  second  integral  Ep, 
the  "error  of  commission".  He  then  found  an  expression  for  Ep,  to  be  used 
for  estimation  of  the  error  of  the  mean  global  error.  In  summary  we  have: 


and 

tit)  ■  — r -  (ID 

s1nJ(*/2) 

Now  we  can  expand  N  in  a  series  of  surface  harmonics  of  degree  n,  as: 


N  ■  JNn  =  U,i)  =IZN"  U,o) 

n=2  n=2  n=2 


and 


where 


(12) 


(13) 


Ro  ■  4  •WT^o7§}  ■  R,  =  R0  - « [i  -  si"  (W2)]  (14) 

2p>n(g°s  *°)  -  Pp  (cos  Xpo )J 
Rn  =  (n  -  "T )  sin  (*o/2) -  +  2  Rn-1  '  Rn-2 


and  Pn  ( x  )  are  Conventional  Legendre  Polynomials.  Buglia  (1976)  gives 
tables  of  Rn>  recalculated  here  in  Appendix  A. 

Now  equation  10  is  an  identity.  That  is,  if  all  the  N^  were  known, 
we  assume  Np  is  known,  then  Ep  could  be  calculated,  subtracted  from  the 
integral  over  the  cap,  and  Agp  would  be  known  without  any  error  from 
neglecting  the  outer  zone.  Of  course  the  N-j  are  not  completely  known,  or 
we  would  not  need  to  continue  this  exercise.  Even  so,  we  can  calculate  the 
R0Np  part  of  Ep.  In  addition  we  note  that  if  some  of  the  Nn  are  known, 
their  contribution  to  Ep  can  be  computed.  We  could  assume  the  lower 
degree  coefficients  of  the  geopotential  are  known  with  sufficient  accuracy 
to  adopt  them  as  known.  In  this  case  the  simplest  procedure  is  to  sub¬ 
tract  this  known  geopotential  surface  of  degree  and  order  nQ  from  the 


geo  id  height  data  and  we  have: 


Ep  ■  IT-  (Ro  NP  -  E  Rn  Nn)  (’3) 

V  n*n0+l  / 

to  calculate  the  error  in  neglecting  the  outer  zone. 

Finally  one  can  profit  from  the  variation  of  R  and  N  with  increas- 

n  n 

ing  n.  Recall  that  we  can  write: 


N 


P 


Nn  (  a  ■  0,  if>  =  0) 


(15) 


and  note  that 

OO 

]CRn  Nn 

"="o+1 

is  a  linear  combination  of  the  Nn.  For  example,  if  we  assume  that  present 
geopotential  models  are  sufficiently  accurate  to  degree  10,  and  we  elect 
to  use  a  spherical  cap  with  ’lo=5°,  then  the  relevant  values  of  Rn, 
n=ll,  12,  .  .  .  decrease  slowly  from  22.26617  to  0.029020  at  degree  18. 
However,  we  don't,  in  general,  know  Np  for  each  computation  point. 

Now,  in  general,  if  we  assume  a  geopotential  model  is  sufficiently 
well  known  through  degree  and  order  nQ,  i.e.,  the  unknown  harmonics  are  of 
degree  n  +1  and  higher,  we  can  write: 


where  we  can  compute  the  first  part  of  16,  and  the  remaining  error  Ep  is 
smaller  than  Ep.  Therefore  the  actual  Molodensky  integral  should  be  written 


y  1  -  — 2 N  - 


R  - 

O  nr 


1 1>0 


y 

p  "  16» R 


r/  f 

ip  =o  a  =  o 


f(f)  (N-Np)  sin  'I'd  ^d« 


which  we  call  the  Modified  Molodensky  Integral. 

We  can  now  establish  the  global  mean  square  error  for  Ep  as: 


E(R"- V2  "Kl 


"=V2 

-5  2 

Using  the  celebrated  10  /n  rule  for  the  anomalous  potential  we  have: 


M  jNn2|  =  R2  x  10'10  (2n+l )/n4 


72  x  10"10 


V  ‘Rn  -  \  >2  122J1L 


(21) 


Table  2  shows  representative  values  for  this  error  term. 

The  case  nQ=0  is  adopting  no  reference  geopotential  model.  For  prac¬ 
tical  limitations  of  a  spherical  cap  radius  of  20°  or  less,  the  error  of 
omitting  the  outer  zone  exceeds  several  hundred  mill igal s . 

It  is  customary  to  assume  that  the  reference  ellipsoid  is  correctly 
positioned  with  respect  to  the  earth's  center  of  mass  and  principal  moments 
of  inertia.  In  this  case,  one  actually  assumes  knowledge  of  the  reference 
surface  nQ=2,  and  the  Modified  Molodensky  Integral  provides  an  accuracy  of 
several  tens  of  milligals.  In  fact  there  would  be  terms  of  degree  and  order 
2  remaining,  and  a  more  careful  calculation  should  be  made  for  the  case  of 
adopting  a  best  fitting  ellipsoid  in  contrast  to  a  reference  surface  of 
degree  two. 

However,  adopting  a  reference  geopotential  of  degree  and  order  12,  and 
using  a  spherical  cap  of  6°  would  give  a  1.24  mill igal  gravity  anomaly. 

Table  2  can  be  used  as  a  prescription  for  choosing  the  cap  size  and 
reference  geopotential  model  for  any  desired  gravity  anomaly  accuracy,  due 
to  neglecting  the  outer  zone  in  the  integral.  It  also  dramatically  illus¬ 
trates  the  benefits  of  using  an  accurate  low  degree  and  order  reference  geo¬ 
potential  model.  With  an  accurate  geopotential  model  of  degree  and  order 
18  and  a  six  degree  integration  cap,  the  error  due  to  ignoring  the  outer 
zone  would  be  reduced  to  0.07  milligals. 


MODIFIED  MOLODENSKY  FORMULA  ERROR  (MGALS) 


IV. 


AN  ALGORITHM  FOR  THE  MOLODENSKY  INTEGRAL 


The  practical  implementation  of  the  Molodensky  Integral  investigated 
here  is  taken  from  a  computer  code  developed  at  DMAAC  (St.  Louis)  and  kindly 
provided  to  MATHEMATICAL  GEOSCIENCES  INC.  for  this  study. 

The  Molodensky  program  assumes  a  data  base  exists,  with  complete 
coverage  of  the  region  of  interest.  The  spherical  region  is  gridded  in 
latitude  and  longitude,  resulting  in  rectangular  areas.  The  gridding  is 
uniform,  and  for  this  purpose  was  chosen  to  be  15  minute  by  15  minute  rec¬ 
tangles.  The  region  selected  for  study  was  at  approximately  60°  latitude. 
Therefore  the  areas  were  approximately  28km  x  14km  rectangles. 

The  Molodensky  program,  takes  as  input  the  region  containing  the  com¬ 
putation  points,  the  spacing  of  the  computation  points,  and  the  cap  size. 
Then,  for  each  computation  point,  all  the  points  in  the  spherical  cap  are 
extracted  from  the  data  base  and  the  quadrature  is  performed.  The  rectangle 
containing  the  computation  point  (called  the  inner  zone)  is  treated 
separately,  necessitated  by  the  singularity  of  the  Kernal  (11)  at  the  com¬ 
putation  point.  This  square  is  subdivided  with  a  uniform  6x6  grid.  A  two 
dimensional  bi-cubic  spline  is  fit  to  the  25  points  closes  to  the  computa¬ 
tion  point.  The  spline  fit  is  used  to  calculate  the  geoid  values  on  the 
6x6  grid.  These  36  sub  areas  are  then  integrated  in  the  same  way  as  the 
other  rectangles.  The  treatment  of  the  inner  zone  will  be  discussed 
further  in  Section  V.  The  actual  integration  takes  the  value  of  the  Kernal 
at  the  center  of  the  rectangle,  and  multiplies  by  the  area  of  the  region. 

As  will  be  seen,  this  quantitization  works  satisfactorily. 
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This  program  has  been  modified  in  a  number  of  ways.  First,  for  test¬ 
ing  purposes,  the  data  base  was  expanded  to  include  a  computed  gravity 
anomaly  at  each  grid  point.  When  developing  simulated  data  for  testing, 
the  gravity  anomaly  was  also  calculated  and  included  in  the  data  base.  The 
program  was  modified  to  retrievethe  gravity  anomaly  and  compare  the  value 
computed  from  the  Molodensky  Integral  with  the  simulated  value.  Next,  the 
program  was  modified  to  calculate  and  explicitly  use  the  error  function 
(Rn,  eq  14)  and  explicitly  use  the  Modified  Molodensky  Integral  (18). 
Finally,  alternate  methods  to  treat  the  inner  zone  were  developed  (see 
Section  V). 

The  simulated  data  base  was  developed  using  potential  coefficients. 

In  addition,  noise  was  added  to  the  data.  The  first  tests  were  done  to 
show  the  whole  calculation  was  self-consistent.  For  this,  recall  that  the 
Modified  Molodensky  Integral  has  an  error  term  (17)  that  involved  surface 
harmonics  of  degrees  greater  than  the  reference  field  employed.  For  this 
simulation,  only  harmonics  of  degree  twelve  were  used.  Therefore,  the  error 
term,  for  higher  harmonics,  should  be  exactly  zero.  The  root  mean  square 
agreement  of  gravity  anomalies  estimated  with  the  Modified  Molodensky 
formula,  and  the  original  values  was  0.047  mgals.  This  departure  from  zero 
is  consistent  with  the  accuracy  of  the  geoid  values  and  gravity  anomalies 
stored  in  the  data  base  (1  millimeter,  and  10  microgals).  The  same  data 
was  corrupted  with  10  cm  Gaussian  Noise.  In  this  case  the  standard  error 
was  14.958  milligals.  The  results  are  summarized  in  Table  3. 
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table  3 

MODIFIED  MOLODENSKY  PROGRAM  ACCURACY 
REFERENCE  FIELD  12TH  DEGREE  ONLY 


sigma  on  data 
(cm) 

<0.1 

10.0 


sigma  (Ag) 
(mgal ) 

0.047 

14.958 


The  fact  that  the  zero  noise  case  agreed  so  well  demonstrated  that  the  basic 
theory  and  implementation  work.  However,  the  significant  disagreement  of 
the  estimate  with  Table  1  mandates  further  investigation.  In  particular 
the  treatment  of  the  inner  zone  and  the  geoid  value  at  the  computation  point 
have  been  found  to  be  critical.  The  original  code  performed  a  spline  fit  to 
the  nearest  25  points  to  the  computation  point.  It  determined  64  spline  co¬ 
efficients,  i.e. ,  more  unknowns  than  observations.  Therefore  the  fit  was 
essentially  an  interpolation  exactly  matching  the  observed  values,  including 
the  computation  point.  Any  error  at  the  computation  point  directly  corrupts 
every  point  in  the  integration  since  the  integration  involves  N-Np.  We 
therefore  explored  the  treatment  of  the  inner  zone,  to  improve  the  value 
N  and  to  improve  the  treatment  of  the  inner  zone  integral. 
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V. 


COMPUTATION  OF  THE  INNER  ZONE 

V.l .  Use  of  the  Bicubic  Spline  Fit 

As  discussed  above,  the  spline  fit  was  made  to  interpolate  the  25 
points  closest  to  the  computation  point.  Its  primary  use  is  for  calculation 
of  N  near  to  Np  on  a  grid  finer  than  the  data  are  given.  The  interpolation, 
in  fact,  must  go  through  the  computation  point  because  the  kernal  is 
singular  at  the  computation  point.  Put  another  way,  near  the  computation 
point  N-Np  must  approach  zero  faster  than  sin  (<fc/2)  or  the  integration  is 
singular. 

The  number  of  data  points  used  for  the  spline  fit  was  extended  from  25 
to  81.  The  calculation  then  had  fewer  unknowns  than  data  points,  which  re¬ 
sults  in  a  kind  of  Least  Squares  fit  to  the  data.  However,  the  spline  still 
fit  exactly  at  the  computation  point,  and  did  not  allow  an  improved  value 
of  Np.  However,  the  values  of  N  used  in  the  inner  zone  improved  the  fit  to 
12.198  milligals. 

To  further  quantify  the  effect  of  errors  in  Np  and  the  inner  zone,  on 
the  gravity  anomaly,  the  true  values  of  N  were  used,  i.e.,  noise  was  not 
added  to  the  values  used  for  calculating  the  spline  fit  and  interpolation  in 
the  inner  zone.  Therefore  in  the  integration,  the  true  value  of  Np  and  the 
inner  zone  integral  were  used.  In  this  case  the  fit  was  improved  to  2.884 
milligals.  Of  course,  this  could  never  be  achieved  in  practice,  since  one 
never  knows  the  true  value  of  Np.  This  shows,  however,  the  absolute  neces¬ 
sity  for  optimum  treatment  of  the  geoid  values  at  the  computation  point. 
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This  demonstration  that  the  inner  zone  treatment  is  a  primary  error  source 
in  the  quadrature  led  to  alternate  methods  of  treating  the  inner  zone. 

V.2.  Least  Squares  Fit  for  the  Inner  Zone 

The  weakness  of  the  spline  fit  is  that  it  does  not  provide  a  Least 

Squares  estimate  of  the  polynomial  coefficients,  Np,  and  N  within  the  inner 

zone.  We  therefore  sought  a  fit  to  a  number  of  values  of  N  near  to  N  .  We 

P 

selected  the  simple  relation: 

J  —  3  j  j  2  2 

N  (x ,  y)  =  ^  x  y  +  a£2  x  y  (23) 

i-j  J=o 

where  x,  y  are  coordinates  with  respect  to  the  computation  point.  There¬ 
fore: 


Np  =  N(o,o) 


oo 


(24) 


This  would  provide  the  desired  Least  Squares  estimate  of  Np.  As  a  second 
benefit,  we  can  directly  estimate  the  inner  zone  integral  from  the  coeffi¬ 
cients  a... 

*  J 

If  we  assume  d<r  is  small,  then  in  the  limiting  flat  earth  approxima¬ 
tion,  the  Molodensky  Kernal  becomes: 


f  (4) 


1 


sin  (<|>/2) 


2  °  2  3/? 

(*  +  y2)  2 


(25) 
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Now,  consider  the  integral  in  (18): 


4o  2*  40  2 * 


4=0  ot- o  4=0  a=o 


(26) 

The  only  non  zero  integrals  are: 

A2,0  =  A0,2  =  *  »  A22  */r|2  (27) 

So  the  integral  for  the  inner  zone  is: 

Ag  inner  =  "  ~ ?R~  j^a20  +  a20)  ro  "  a22  r03/12J 

This  relation  illustrated  the  nature  of  the  Molodensky  integral. 

The  Stokes  Integral  (1)  to  obtain  the  geoid  (N)  from  gravity  anomalies 
(Ag)  is  essentially  a  smoothing  operation.  One  can  tell  this  naively  from 
the  fact  that  a  geoid  map  is  much  smoother  than  a  gravity  map.  The  inverse 
process  to  obtain  gravity  (ag)  from  the  geoid  (N)  must  be  the  inverse,  or 
a  roughing  operation.  Equation  (29)  shows  that  the  Molodensky  Integral  es¬ 
sentially  uses  the  second  degree  curvature  of  the  geoid  to  infer  gravity, 
i.e.,  2nd  and  higher  derivatives.  It  also  points  out  what  we  already  know, 
that  the  absolute  value  of  the  geoid  doesn't  play  a  role  in  estimating 
gravity,  only  the  derivatives.  The  roughing  character  of  the  Molodensky 
Integral  must  therefore  amplify  short  wavelength  variations  in  N,  which 
is  serious  if  these  short  wavelength  variations  are  errors. 


With  this  Least  Squares  fit  for  the  inner  zone,  a  number  of  experi¬ 
ments  can  be  performed.  First,  the  fit  was  used  with  the  noiseless  data. 

In  this  case,  the  fit  was  0.015  milligals  ms.  This  improvement  over  the 
spline  fit  interpolation,  0.047  mgals,  was  due  to  improvement  of  Np  at  the 
submillimeter  level. 

To  test  the  analytical  integral  of  the  inner  zone  (28),  one  can  do  two 
things.  The  inner  zone  is  a  rectangle,  in  fact  it  is  a  trapezoid  but  far 
from  the  pole  this  error  is  small  and  is  ignored,  and  (28)  is  for  a  circle. 
One  can  choose  the  circle's  radius  such  that  it  has  the  same  area  as  the 
rectangle,  i.e.: 


Vab 
* 

Alternatively  one  can  choose  r0~b,  the  smaller  dimension,  and  integrate  the 
remaining  area  of  the  rectangle  as  before.  In  this  case  one  selects  areas 
in  the  inner  zone  to  include  in  the  integral,  those  areas  with  a  mid  point 
outside  the  radius  r0-  There  is  an  error  associated  with  either  approach. 
In  both  cases  more  careful  calculation  could  be  developed,  but  this  doesn't 
seem  worth  pursuing  given  the  quality  of  the  fit.  In  both  cases  the  fit  to 
the  noiseless  data  was  0.018  milligals,  verifying  the  analytical  integral. 
This  small  difference  between  the  direct  numerical  quadrature  of  the  6x6 
inner  zone  grid  and  the  analytical  formula  is  also  true  for  treating  data 
with  Gaussian  Noise.  Therefore  we  adopted  the  direct  numerical  quadrature 
for  actual  calculation. 

Finally,  the  Least  Squares  polynomial  fit  was  used  to  obtain  gravity 


anomalies  from  data  with  10cm  Gaussian  noise.  The  rms  came  to  5.033 


milligals,  which  is  in  satisfactory  agreement  with  4.55  mi 11 i gal s  from 
Table  1.  The  results  of  all  these  numerical  experiments  are  summarized  in 
Table  4. 


TABLE  4 

ACCURACY  OF  MODIFIED  MOLODENSKY  INTEGRAL 
REFERENCE  FIELD  DEGREE  12  ONLY 


Inner  Zone 

Noise  on  Data 
(cm) 

Sigma  (Ag) 
(mqal ) 

Spline  to  25  points 

<0.1 

0.047 

Spline  to  25  points 

10.0 

14.958 

Spline  to  81  points 

10.0 

12.198 

Spline  to  81  points 
(no  inner  zone  noise) 

10.0 

2.884 

Least  Squares  Polynomial 

<0.1 

0.015 

Least  Squares  Polynomial 
(Analytical  Inner  Zone) 

<0.1 

0.018 

Least  Squares  Polynomial 

10.0 

5.033 

VI. 

SUMMARY  OF  MOLODENSKY  INTEGRAL 

The  numerical  calculation  of  gravity  anomalies  from  geoid  data  (as  ob¬ 
tained  for  example,  from  satellite  to  sea  surface  altimetry)  has  been  dis¬ 
cussed  in  terms  of  the  Molodensky  Integral.  Four  sources  of  error  have 
been  investigated: 

1)  Error  from  the  neglected  outer  zone 

2)  Error  from  singularity  of  Molodensky  Kernal 

3)  Error  from  Quantitization  due  to  replacing  an  integral  with 
a  summation 

4)  Error  from  Gaussian  noise  on  the  data. 

It  is  found  that  error  2)  and  3)  are  negligible. 

The  use  of  a  spherical  cap  is  imposed  in  two  ways.  First, as  a  practical 
matter  performing  a  quadrature  over  the  whole  unit  sphere  is  prohibitive, 
and  second,  with  a  satellite  to  sea  surface  altimeter,  one  does  not  have 
geoid  data  everywhere.  A  Modified  Molodensky  Integral  (18)  is  given  that 
reduces  the  error  in  neglecting  the  outer  zone.  Further  significant  improve¬ 
ments  can  be  made  by  using  a  reference  geopotential  model  for  the  low  degree 
and  order  harmonics.  Table  3  gives  the  expected  errors  in  gravity  anomalies 
for  a  given  spherical  cap  size  and  reference  geopotential  models.  Today 
geopotential  models  are  quite  accurate  to  degree  and  order  10,  and  signifi¬ 
cant  improvements  may  be  available  in  the  next  five  years.  The  benefits  in 
using  a  reference  geopotential  model  are  enormous. 

Numerical  experiments  using  simulated  noisy  data,  confirm  the  a  priori 

accuracy  assessments.  Table  1.  For  example,  with  10cm  geoid  values  a 

1 5 ' xl 5 '  gravity  anomaly  can  be  obtained  with  an  accuracy  of  approximately 
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5  milligals.  To  obtain  this  accuracy  it  is  essential  to  treat  the  geoid  at 
the  computation  point  carefully.  This  leads  to  the  second  conclusion  that 
special  care  must  be  taken  to  achieve  the  most  accurate  geoid  values.  In 
particular  to  reach  one  milligal  for  1 5 ' xl 5 1  areas  one  must  have  two  cm 
geoid  values. 

The  implications  for  processing  the  Seasat-1  and  Geosat-A  data  are 
clear.  Great  care  must  be  taken  for  all  error  sources.  The  long  wavelength 
orbit  errors  are  not  critical  as  long  as  the  orbit  bias  has  been  removed  with 
a  crossing  arc  analysis  or  an  orbit  matching  adjustment  such  as  the  spline 
fit  used  by  Cloutier  et  al  (1981).  However,  any  such  fit  must  be  achieved 
with  centimeter  accuracy.  In  addition,  parasitic  errors  such  as  refraction, 
real  oceanographic  signals,  and  instrument  noise  must  be  eliminated  to  the 
extent  possible.  The  lack  of  a  radiometer  on  Geosat-A  such  as  was  available 
on  Seasat-1  could  create  a  serious  lack  of  important  data.  Approaches  to 
these  problems  will  be  discussed  in  the  next  part  of  this  Report. 


INSTRUMENT  AND  ORBIT  ERROR 


The  accuracy  of  the  altimeter  is  fixed.  The  Geosat-A  altimeter  is  in 
tended  to  be  identical  to  that  on  Seasat-1.  The  characteristics  are  sum¬ 
marized  in  Table  5.  There  are  several  instrument  related  corrections 
necessary.  For  example,  the  significant  wave  height  or  sea  state  bias  has 
been  analyzed  (Born,  Richards,  &  Rosborough,  1982;  Hayne  &  Hancock,  1983; 
West,  1981).  None  of  these  analyses  give  formal  statistics.  However,  it 
is  estimated  that  the  uncertainty  in  the  sea  state  bias  is  2%  of  the  signi 
ficant  wave  weight.  Therefore  for  hl/3  <  2.5  m,  one  would  have  less  than 
5cm  rms  error.  There  is  also  the  response  of  the  phase  lock  loop  in  the 
altimeter  tracker.  In  regions  of  rapid  change  of  altitude,  this  is  an  im¬ 
portant  correction  (Tapley,  Born,  &  Parke,  1982). 


TABLE  5 

ALTIMETER  CHARACTERISTICS 


Observable 

Altitude  (Seasat) 
Altitude  (Geosat*) 
Significant  Wave  Height 
Swath  Width 


Specification 

8  cm  (precision)  hl/3  <  5m 
3.5  cm  (precision) 

10%  or  0.5m 

2.4  to  12  km  depends  on  seastate 


♦Target  for  Geosat-A,  CDR  November  1982. 


The  absolute  accuracy  of  the  Geosat-A  emphemeris  is  beyond  the  scope 
of  this  investigation.  There  are  a  number  of  approaches  that  can  be  used 
to  optimize  it.  Of  course,  some  improvement  of  the  orbit  accuracy  over  that 
for  Seasat-1  will  naturally  come  from  two  things:  1)  Geosat-A  is  physically 
more  compact,  without  the  large  antennae  that  were  on  Seasat-1.  This  will 
reduce  the  errors  in  modelling  solar  radiation  pressure  and  atmospheric  drag. 
2)  There  are  improved  gravity  field  models  available,  some  derived  using 
Seasat-1  data,  that  will  also  result  in  more  accurate  ephemerides. 

In  addition,  methods  have  been  developed  to  use  the  observations  at  al¬ 
timeter  track  crossing  points  to  obtain  improved  estimates  of  the  long  wave¬ 
length  orbit  errors.  These  all  use  the  (erroneous)  idea  that  the  sea  sur¬ 
face  height  at  each  crossing  point  is  a  constant.  Therefore,  any  difference 
in  the  observed  height  is  assumed  to  be  an  error  in  the  radial  component  of 
the  satellite  emphemeris.  This  would  be  true  if  there  were  no  time  variable 
oceanographic,  meteorological,  and  observation  biases.  Within  this  assump¬ 
tion,  crossing  arc  analyses  have  achieved  (<10cm),  Cheney  &  Marsh  from  1.5 
years  of  Geos-3  data,  (15  to  22  cm)  Cloutier,  and  (28  cm)  Rapp  with  Seasat-1 
data.  One  would  suppose  these  to  be  optimistic  estimates  since  a  number  of 
the  time  variable  biases  have  periods  longer  than  the  span  of  Seasat-1  data, 

3  months.  Conversely,  one  would  not  expect  the  oceanographic  and  meteorologi 
cal  errors  to  be  correlated  with  orbit  errors.  However,  one  expects  orbit 
errors  to  be  correlated  with  geoid  errors  (Anderle  &  Hoskin,  1977). 

In  summary,  considerable  progress  is  necessary  to  reduce  the  Geosat-A 
orbit  error  to  the  accuracy  of  the  Geosat-A  altimeter.  For  the  initial 
analysis  of  Geosat-A  data,  this  error  source,  after  the  crossing  arc 
analysis  has  been  applied,  should  be  significantly  less  than  20cm.  With 


THE  OCEANOGRAPHIC  FIELD 


VIII. 1  The  Mean  Circulation 

The  principal  non-geodetic  signal  in  the  altimeter  data  comes  from  the 
ocean  circulation  and  eddy  field.  The  mean  ocean  circulation,  if  it  can  be 
defined,  would  be  a  position  dependent,  time  invariant  signal,  and  indis¬ 
tinguishable  from  the  geoid  with  an  altimeter.  It  has  an  amplitude  of  1m 
and  a  length  scale  ranging  from  100  km  to  10000  km.  The  sharpest  gradients 
occur  at  the  strong  western  boundary  currents  such  as  the  Gulf  Stream  or  the 
Kurshio.  At  a  western  boundary  current  a  one  meter  head  rises  in  roughly 
100km.  The  return  flow  is  distributed  across  the  remaining  ocean  basin, 
i.e.,  the  return  flow  represents  a  1  meter  drop  over  10000  km.  Therefore 
the  mean  circulation  is  only  a  factor  near  the  strong  western  boundary  cur¬ 
rents.  Approaches  to  simultaneously  determine  the  geoid  and  the  ocean  cir¬ 
culation  using  altimeter,  gravity,  and  hydrographic  data  have  been  discussed 
(Wunsch,  1978;  Wunsch  &  Gaposchkin,  1980).  At  present  the  very  long  wave¬ 
length  features  of  the  circulation  are  being  investigated  with  Seasat-1 
altimeter  data  using  a  satellite  derived  geoid  as  a  reference  (Tai  &  Wunsch, 
1983).  However,  this  remains  probably  the  most  challenging  and  complex 
interdisciplinary  problem  in  Altimetric  Oceanography. 

VIII. 2  The  Eddy  Field 

Superimposed  on  the  mean  circulation  is  a  variability  with  one  component 
labeled  the  Mesoscale  Eddy  Field.  Many  discussions  make  a  distinction 
between  cold  cyclones  (rings)  and  eddies.  The  former  are  analagous  to 


weather  cyclones  and  are  believed  to  originate  from  similar  dynamical 
phenomena.  There  is  not  yet  any  understanding  of  the  generation  of  eddies. 
They  may  be  rings  that  have  traveled  great  distances,  or  be  generated  by 
unstable  interactions  of  the  flow  and  bottom  topography.  The  flow  field 
associated  with  surface  eddies  can  extend  below  3km.  Rings  are  thought  to 
have  a  lifetime  of  two  to  three  years.  It  is  thought  that  about  five  rings 
are  formed  each  year  from  the  Gulf  Stream.  Therefore  we  expect  to  find 
about  15  rings  in  the  Western  North  Atlantic  at  any  time.  As  emphasized  by 
Koshlyakov  &  Monin  (1978),  the  North  Pacific  is  even  more  complex. 

Rings  and  eddy  fields  are  broadly  characterized  with  a  1  to  2m  ampli¬ 
tude,  a  10  to  30cm  variance,  a  150km  to  250km  horizontal  wavelength,  and  a 
time  constant  of  2  to  5  months  (Fu,  1982;  Koshlyakov  &  Monin,  1978;  Douglas 
&  Cheney,  1981).  By  this  time  constant  we  mean  that  the  variability  will 
move  by  one  wavelength  in  one  time  constant.  Put  another  way,  one  can  hope 
to  average  the  mesoscale  eddy  field  over  one  time  constant.  Unfortunately, 
the  ocean  variability  is  thought  to  have  components  with  decade  time  con¬ 
stants,  and  possibly  longer.  The  oceanographic  data  are  not  yet  sufficient 
to  constain  the  period  of  oceanographic  variability  (Wunsch,  1983).  There¬ 
fore  the  analysis  of  altimetric  data  must  be  hardened  against  this  possible 
low  frequency  oceanographic  variability. 

To  quantify  eddy  field  averaging,  we  can  use  the  description  of  the 
eddy  field  in  terms  of  Rossby  Waves  resulting  from  Baroclinic  Instability. 
The  theory  predicts  the  phase  velocity  of  propagation  c  as: 


0  =  2  Q  cos  <t>  /a£ 


x  is  the  wavelength  (peak  to  trough) 
a  is  the  rotation  rate  of  the  earth  (0.727x10"^  sec-1) 
<Ms  the  latitude 
ag  is  the  earth  radius 


The  baroclinic  Rossby  Radius  of  deformation: 


where 


f  =  2  Q  sin  <t> 


and  g  is  the  acceleration  of  gravity,  h  is  the  height  (head)  of  the  rotating 
system  free  surface.  The  Rossby  Radius  is  the  shortest  distance  over  which 
geostrophic  flows  of  any  significance  can  occur.  At  mid  latitude  for  a 
height  of  1m,  =  31km.  Of  course,  for  a  given  height,  longer  wavelengths 


exist. 


We  can  calculate  the  time  for  a  Rossby  Wave  (eddy)  to  travel  one  wave¬ 


length  as: 


e 

a  cosspX 


$ 


and  for  the  shortest  wavelength  of  a  Rossby  Radius: 

2  ae 

Pn  = - - -  tan  <0 

yfin 

Now,  for  a  mission  duration,  T,  and  with  uniform  sampling,  each  eddy  will 
move  n=T/P  wavelengths.  Therefore  the  error,  after  averaging,  will  be 
h^7  i.e. 

h 

vT~ 

and  for  a  Rossby  Radius  Wave: 


We  give  in  Table  6  some  sample  values,  for  a  mission  of  T=2  years,  and 
sin  <0=cos  <0=1  A/T. 


TABLE  6 


AVERAGING  OF  EDDY  FIELD 
TWO  YEAR  MISSION 
LATITUDE  =  45  DEG. 


h 

X 

error 

Icml 

(km) 

(cm) 

100 

31.3* 

25 

50 

22.0* 

15 

100 

100 

12 

100 

200 

8 

*Rossby  Radius  corresponding  to  h. 

If  we  average  all  the  altimetric  data,  over  as  long  a  time  base  as 
possible,  we  define  a  mean  surface  in  this  way.  The  availability  of  a  10 
year  time  base  using  the,  admittedly  less  accurate,  Geos-3  data  is  an  advan¬ 
tage.  This  average  surface  would  contain  a  mean  circulation  which  could  be 
removed  with  any  number  of  models,  e.g.,  Tai  &  Wunsch  (1983),  the  Steric 
Anomaly  field  of  Wyrtki  (1974,  1975),  or  Levitus  &  Oort  (1977).  For  the 
moment  any  such  model  must  be  viewed  with  caution  as  it  will  inevitably  be 
built  on  some  assumptions.  One  such  assumption  is  the  level  of  no  motion 
which  is  certainly  not  correct;  see  e.g.,  Wunsch  (1978).  From  the  geodetic 
viewpoint,  this  oceanographic  signal  will  be  considered  as  part  of  the 
error  budget  for  some  time  to  come.  Only  continued  fundamental  oceanographic 
research,  an  oceanographic  research  satellite  such  as  TOPEX  (Wunsch,  1981), 
or  a  dedicated  Gravity  Satellite  such  as  GRM  (Kaula,  1983)  will  signifi¬ 
cantly  change  this  situation. 


ENVIRONMENTAL  CORRECTIONS,  A  SURVEY 


The  environmental  contributions  to  the  observations  include: 

hw,  the  wet  tropospheric  correction 

hd,  the  dry  tropospheric  correction 
hb,  the  inverse  barometer  correction 
hi,  the  iohosophere  correction 

he,  the  solid  earth  tide  correction 
ho,  the  ocean  tide  correction 

hs,  the  significant  wave  height  correction 

There  is  general  agreement  that  tide  models  for  he  have  an  error  less 
than  2  cm.  Ocean  tide  models  (ho),  (e.g.,  Schwiderski,  1982a,  b,  c;  Parke 
&  Henderschott,  1980;  Parke,  1980)  are  generally  better  than  10  cm,  where 
applicable.  Parke  &  Henderschott  (1980)  is  primarily  a  deep  sea  model, 
whereas  Schwiderski  (1982a,b,c)  would  be  more  appropriate  near  coasts  where 
tide  gauges  exist.  The  model  accuracy  in  regions  without  tide  gage  measure¬ 
ments  is  unknown.  There  is  therefore  a  possible  true  variable  error  source 
with  a  wavelength  between  500km  and  1000km  (Tapley,  et  al  1983). 

With  suitable  care,  the  ionosphere  correction  can  be  modeled  with  an 
accuracy  better  than  2  cm,  as  described  by  Lorell  et  al  (1983).  Of  course 
a  dual  frequency  altimeter  would  have  measured  this  correction  directly. 

The  inverse  barometer  correction  (hb)  corresponds  to  the  depression 
(approximately  lcm/mbar)  of  the  sea  surface  with  change  in  atmospheric 
pressure.  This  is  the  order  of  25  cm  with  a  time  constant  of  approxi¬ 
mately  two  days,  i.e.,  the  time  constant  of  oceanospheric  response  is 
about  48  hours.  Except  near  storms  or  other  transient  phenomena,  the 


FNOC  averaged  global  pressure  field  will  provide  the  necessary  correction. 

The  dry  troposphere  correction  (hd)  has  an  amplitude  of  approximately 
2.1  m.  The  surface  pressure  available  from  the  FNOC  is  adequate  to  compute 
this  correction  with  1  or  2  cm  accuracy.  In  fact  a  comparison  of  FNOC  data 
with  radiosonde  data  save  an  agreement  of  0.59  cm  rms  for  this  dry  tropo¬ 
sphere  correction  (Tapley  et  al  1983). 

The  wet  troposphere  correction  (hw)  can  exceed  40  cm,  with  a  length 
scale  of  150  to  500  km.  On  Seasat-1 ,  the  Scanning  Multichannel  Microwave 
Radiometer  (SMMR)  provided  data  to  compute  this  correction  with  an  RMS  of 
2.7  cm  rms.  The  use  of  FNOC  data  for  computing  this  correction  is  the  sub¬ 
ject  of  the  remainder  of  this  report.  The  accuracy  that  can  be  obtained 
using  FNOC,  averaged  over  several  months,  can  be  approximately  2  cm,  and 
will  not  preclude  obtaining  a  mean  ocean  surface  with  an  accuracy  better 
than  5  cm. 

The  accuracy  of  the  SMMR  wet  troposphere  correction  (hw)  is  2.79  cm 
rms,  as  obtained  in  comparison  with  radiosonde  data  (Tapley  et  al  1982b). 
The  FNOC  derived  data  was  also  discussed  and  shown  to  have  an  accuracy  of 
5.73  cm  rms.  The  SMMR  had  a  horizontal  resolution  of  about  40  km.  The 
FNOC  information  is  based  on  hindcast  models  with  a  grid  spacing  of  several 
hundred  kilometers.  Furthermore  data  on  which  the  for  FNOC  models  are 
based  is  less  complete  in  the  southern  hemisphere.  One  would  expect  the 
FNOC  derived  correction  to  be  less  accurate.  The  remainder  of  this  report 
is  concerned  with  quantifying  the  wet  troposphere  correction  further,  and 
examining  how  well  one  can  expect  to  make  this  correction  in  the  absence 
of  SMMR,  or  equivalent,  data. 


TROPOSPHERIC  REFRACTION 


X.l  Basic  Principles 

The  increase  in  time  delay  of  electromagnetic  waves  depends  on  the 
index  of  refraction  n.  The  velocity  (v)  of  propagation,  v= c/n,  where  c  is 
the  speed  of  light.  In  the  case  of  the  troposphere  the  index  of  refraction 
is  nondispersive  for  frequencies  below  the  22.5  GHz  water  vapor  absorption 
line.  The  Seasat  altimeter  operates  at  a  single  frequency  of  13.5  GHz. 
Therefore  the  group,  phase,  and  signal  velocities  are  all  the  same  and  are 
less  than  c. 

The  index  of  refraction  can  be  written  (Bean  &  Dutton,  1966); 

n  =  1  +  N 


where 


N  =  (1-n) 


77.6 


+  3.73  x  10"5 


x  10"6 


where;  z  is  the  height  in  question,  P(z)  is  the  atmospheric  pressure  in 
millibars,  T(z)  is  the  temperature  in  degrees  Kelvin,  e(z)  is  the  partial 
pressure  of  water  vapor  in  millibars,  and  the  constants  77.6  and 
3.73x10**5  depend  on  the  average  global  atmospheric  composition  near  the 
earth's  surface.  The  total  time  delay  of  a  radar  pulse,  assuming  vertical 
propagation,  is: 


AT*/  dz  =  1 
V  C 


n  dz  = 


u_ 0 + 


N)  dz 


Since  the  troposphere  extends  to  altitude  Zt,  the  measured  time  delay  con¬ 
verted  to  distance  h=c  (dt) 


h  = 


/"dz  +  ft 

n  J  n 


where  the  first  term  is  the  geometric  distance,  i.e.,  based  on  the  vacuum 
delay  (hv)  and  the  second  term  (ht)  is  the  additional  delay  due  to  the 
troposphere.  Now: 


1C6  ht  =  77.6 


l  L 

f  -rffj-  dz  +  3-73  *  '°5J  dz 


and  ht=hd+hw.  The  first  term  is  known  as  the  dry  tenrs,  and  is  evaluated 
assuming  the  atmosphere  obeys  the  perfect  gas  law  and  is  in  hydrostatic 
equilibrium,  as: 


2.277  x  10' 


1  +  0.0026  cos 


<2*>]  Po 


where  PQ  is  the  atmospheric  pressure  at  the  surface  and  <p  is  the  latitude 
(geodetic).  The  second  term  accounts  for  the  variation  of  surface  gravity 
Using  the  same  assumptions,  that  the  temperature  decreases  linearly  above 
the  ocean,  and  that; 


IM 

To 
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The  second  term  can  be  integrated  to: 


h  =  2.277x10 
w 


-3 


1255 


0.05 


where  eQ  and  Tq  are  respectively,  the  partial  pressure  water  vapor  in 
millibars  and  temperature  of  the  atmosphere  in  degrees  Kelvin  at  the 
surface. 

The  SMMR  on  Seasat  recorded  antenna  temperatures  at  18,  21,  and  37 
GHz  with  both  horizontal  and  vertical  polarization.  Based  on  the  develop¬ 
ment  of  Wisler  &  Hollinger  (1977)  the  wet  troposphere  range  correction  can 
be  expressed  as: 


6 

R  -  Ro  +  *  £  “i  (Bi  -  B°> 

i=l 

where  R  is  the  range  correction,  are  constants,  are  the  brightness 
temperatures,  B°  are  the  reference  brightness  temperatures ,  RQ  is  the 
reference  range  correction,  and  y  is  an  adjustable  scale  parameter. 

In  the  Seasat-1  case  «.  and  B°  were  adopted  from  Hollinger  &  Wisler 
and  the  (2)  parameters  RQ  and  y  were  determined  using  radiosonde  data  ob¬ 
tained  at  times  coincident  with  Seasat-1  overflights  (Tapley  et  al  1982). 
Sixteen  (16)  values  from  four  geographical  locations  (Kwajalin,  Bermuda, 
and  Ships  Endurer  and  Hecla)  were  used  giving  R0=11.73  (±0.12)  cm  and 
y  =1.29  (±0.0008)  with  a  residual  standard  deviation  of  2.79  cm.  It 
should  be  emphasized  that  the  data  are  small  in  number  (16  points)  and 


that  the  observations  had  either  small  or  large  for  atmospheric  water 
vapour:  there  is  a  lack  of  values  in  the  mid  moisture  range.  Finally, 
the  data  from  a  fifth  site  (Majuro)  were  deleted  because  they  save  incon¬ 
sistent  results.  Since  more  radiosonde  data  corresponding  to  Seasat-1 
tracks  is  now  available,  the  SMMR  calibration  of  water  vapor  could  be  dis 
cussed  more  fully.  If  the  Seasat-1  data  have  anomalies,  apparently  con¬ 
nected  with  wet  troposohere  corrections,  this  would  be  a  fruitful  avenue 
of  investigation.  In  any  event  we  proceed  assuming  that,  where  SMMR  data 
is  available,  the  water  vapor  correction  (hw)  is  known  with  an  accuracy 
of  2.79  cm  rms.  In  the  same  report  (Tapley  et  al  1982),  a  comparison  of 
13  FNOC  derived  hw  corrections  indicate  they  have  an  accuracy  of  5.72  cm 
rms. 

In  what  follows  we  compare  all  SMMR,  h(SMMR),  and  FNOC,  h(FNOC),  cor 
rections  assuming  the  SMMR  values  are  correct.  There  are  three  questions 
in  mind: 

1)  On  average  how  good  is  the  FNOC  correction? 

2)  Are  there  regional  variations? 

3)  How  good  is  the  average  variation  when  used  for  the  correction? 

X.2  The  Seasat  SMMR  Data  Analyzed 

All  available  Seasat  data  in  the  western  North  Atlantic  (20N  -  60N, 
270W  -  330W)  was  requested  from  the  NOAA/EDIS.  This  data  was  provided  on 
one  1600  BPI  magnetic  tape  as  described  in  Lorell  et  al  (1980).  All  the 
h(SMMR)  and  h(FNOC)  data  with  the  latitude  and  longitude  were  culled  from 
the  data.  Only  observations  where  both  values  were  provided  were  used. 
The  resulting  132298  observations  were  used  for  comparison.  The  h(SMMR) 


ard  h(FNOC)  agreed  to  4.097  cm  rms  with  a  mean  difference  of  0.040  cm. 

Thus  to  first  order  one  can  conclude  that  the  h(SMMR)  correction  and 
average  meteorological  data  are  in  excellent  agreement.  Further,  the  use  of 
h(FNOC)  in  the  Western  North  Atlantic  is  quite  acceptable.  If  we  adopt 
2.79  cm  rms  (Tapley  et  al ,  1982)  as  the  accuracy  of  h(SMMR)  and  assume  er¬ 
rors  in  the  two  data  sets  are  uncorrelated,  then  the  h(FNOC)  would  have  an 
accuracy  of  2.94  cm  rms. 

X.3  The  Regional  Variation  of  Wet  Troposphere  Correction 

For  the  use  of  satellite  to  sea  surface  altimetry  for  determination  of 
the  mean  sea  surface  will  use  all  available  data.  Therefore,  one  is  in¬ 
terested  in  the  geographical  variation  of  this  correction,  its  error,  and 
the  autocorrelation  function  of  the  error.  In  other  words,  we  ask  what  the 
error  in  the  average  sea  surface  height  would  be  using  h(FNOC),  assuming 
the  h(SMMR)  to  be  correct,  as  a  function  of  position  and  what  is  its  auto¬ 
correlation  function.  To  address  this  question,  a  data  base  was  formed 
taking  the  average  SMMR  correction  and  the  average  FNOC  correction  for  a 
grid.  A  15'  x  15'  (27km  x  27km)  grid  was  selected.  This  grid  is  more 
coarse  than  the  altimeter  data  spacing  (7  to  15  km),  finer  than  the  SMMR 
resolution  (40  km)  or  the  FNOC  resolution  (250  k.n),  and  approximately  the 
desired  resolution  of  the  final  geodetic  product.  This  data  base  was  used 
for  the  remainder  of  this  analysis. 

To  discover  any  regional  variation  in  hw,  the  mean  SMMR  data  was  again 
averaged  in  2deg  x  2deg  squares.  This  was  a  weighted  average  based  on  all 
observations  in  the  2deg  x  2deg  region  rather  than  a  simple  mean  of  the 
15'  x  15'  values.  Figure  lisa  number  map  of  these  values.  Also  on  the 


number  map  is  a  rough  outline  of  the  North  American  Continent.  There  are 
altimeter  measurements  from  The  Great  Lakes,  St.  Lawrence  River,  Hudson's 
Bay  and  other  water  surfaces.  These  observations  are  small  in  number  and 
do  not  affect  the  overall  statistics.  The  agreement  of  the  FNOC  and  SMMR 
data  for  these  overland  regions  is  good.  Of  course,  these  values  will  not 
contribute  to  determination  of  the  marine  geoid  because  these  bodies  of 
water  are  not  at  the  same  potential  as  mean  sea  level.  Evidently  there  is 
a  broad  latitude  dependence  in  the  wet  troposphere  correction,  with  a  less 
pronounced  longitude  dependence.  The  latitude  dependence  accounts  for  94% 
of  the  variance  in  the  h(SMMR).  For  each  15'  band  of  latitude,  the  mean 
was  taken.  Figure  2  is  a  plot  of  the  latitude  dependence.  The  latitude 
dependence  of  the  h(FNOC)  is  virtually  identical.  Furthermore,  for  each 
latitude  the  variance  of  the  difference  between  the  signal  and  the  mean  was 
calculated  and  is  plotted  in  Figure  3.  This  variance  represents  the  wet 
troposphere  information  not  included  in  the  mean  given  in  Figure  2.  The 
agreement  is  between  3  and  6  cm  rms,  indicating  that  using  a  latitude  de¬ 
pendent  mean  hw  would  give  troposphere  corrections  with  an  accuracy  of  5  cm 
rms.  Figure  3  is  plotted  on  the  same  scale  as  Figure  2  to  demonstrate  the 
remarkable  fit  of  the  mean  troposphere  correction.  Also,  we  see  that  the 
FNOC  data  is  sufficiently  good  to  determine  this  latitude  dependent  func¬ 
tion,  at  least  in  the  Western  North  Atlantic. 


Figure  1.  Number  Map  of  Wet  Troposphere  Correction,  Western  North 
Atlantic  —  2  deg  x  2  deg  means  (cm) 


Figure  2.  Mean  Wet  Troposphere  Correction,  15'  Latitude  Bands 
Western  North  Atlantic 
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X.4  Spectra  of  Wet  Troposphere  Correction 

Finally,  we  consider  the  variance  and  autocovariance  of  the  SMMR  and 
FNOC  derived  tropospheric  corrections  and  their  difference  with  respect  to 
the  mean  latitude  dependent  hw  correction  (H=H(  <*>)).  The  h(SMMR)  variance 
represents  the  true  information  hw  and  would  be  an  estimate  of  the  error 
committed  ignoring  the  hw  correction.  The  h(SMMR)-H  variance  represents 
the  error  in  H.  The  n(FNUC)-H  represents  tne  amount  of  short  wave¬ 
length  information  actually  present  in  the  FNOC  data.  Finally,  the  variance 
of  the  difference,  h(FNOC)-h(SMMR) ,  would  represent  the  error  in  mean  sea 
surface  if  the  FNOC  correction  were  to  be  used.  The  basic  results  are  sum¬ 
marized  in  Table  7. 


TABLE  7 

STATISTICS  OF  SMMR  AND  FNOC  DATA 


<h(SMMR)**2> 
<h(FNWC)**2> 
<(h(SMMR)-H)**2> 
<(h(FN0C)-H)**2> 
<(h(FNOC)-h(SMMR) )**2> 


506.  cm**2 

491 .  cm**2 

20.03  cm**2 
5.29  cm**2 
15.53  cm**2 


(22.5cm)**2 
(22.2cm)**2 
(4.48cm)**2 
(2.3  cm)**2 
(3.94cm)**2 


The  autocovariance  functions  are  presented  in  Figures  4,  5,  and  6.  The 
tabular  data  are  given  in  Table  8.  The  small  value  of  < (h(SMMR)-h(FNOC) ) 
**2>  is  consistent  with  a  meteorological  hindcast  with  coarse  grid  spacing. 
It  also  reemphasizes  the  notion  that  mean  regional  values  are  almost  as 
good  as  the  detailed  FNOC  or  SMMR  values  in  a  statistical  sense.  The 
very  rapid  decrease  in  the  h(SMMR)  autocorrelation  function  indicates  that 
when  spatially  averaged,  the  departure  from  the  mean  (H)  is  very  close  to 
a  random,  uncorrelated  signal.  The  error  in  the  FNOC  is  only  slightly 
less  than  this  ignored  signal. 


h(SMMR) 


Autocorrelation  Function 
Variance  2  0.03  cm2 


Figure  4.  Autocorrelation  Function  for  h(SMMR) 
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Figure  5.  Autocorrelation  Function  for  h(FNOC) 
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Figure  6.  Autocorrelation  Function  for  h(FNOC)-h(SMMR) 


These  results  are  an  average  for  all  15*  x  15'  squares  containing  one 
or  more  data  points.  This  data  set  was  obtained  in  approximately  90  days. 
Some  squares  have  many  more  observations  than  others,  and  we  can  inquire 
if  the  uncertainty  improves  as  the  number  of  data  points  increases.  Figure 
7  plots  the  rms  error  as  a  function  of  the  number  of  data  points  in  the 
mean.  By  far  the  largest  number  of  squares  had  1 ,  2  or  3  data  points,  and 
these  dominate  the  overall  sigma.  However,  for  squares  with  15  or  more 
data  points,  the  sigma  is  approximately  (2.25  cm)**2.  There  does  not  seem 
to  be  much  further  improvement  as  the  number  of  data  points  increases.  It 
is  unclear  why.  However,  with  the  two  year  lifetime  of  Geosat-A  one  can 
reasonably  expect  to  have  all  squares  observed  more  than  15  times.  Of 
course,  the  actual  distribution  of  ground  tracks  depends  on  the  details  of 
the  satellite  manoeuvers  and  on  other  programmatic  considerations.  From  the 
point  of  view  of  averaging  the  wet  troposphere  error  one  should  obtain  as 
uniform  a  geographical  distribution  as  possible. 
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Figure  7.  Error  of  15'  x  15'  Square  as  Function  of  Number  of  Observations 
in  Average. 
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X.5  Discussion  of  Wet  Troposphere  Error 


We  have  taken  the  Seasat-1  data  set  in  the  Western  North  Atlantic  to 
compare  the  wet  troposphere  correction  derived  from  the  SMMR  and  that 
derived  from  the  FNOC  data.  The  SMMR  data  is  believed  to  have  an  accuracy 
of  2.67  cm  rms  and  is  used  as  the  standard  of  comparison.  The  objective 
here  is  to  project  the  situation  for  Geosat-A  which  will  not  be  equipped 
with  a  SMMR,  and  will  therefore  depend  on  FNOC  for  this  data. 

This  study  is  limited  in  three  ways.  First,  it  assumes  the  Seasat-1 
SMMR  is  a  suitable  standard  of  comparison.  Recall  that  the  SMMR  calibra¬ 
tion  was  made  with  only  16  radiosonde  data  sets,  16%  of  the  data  was  rejected, 
and  even  so  had  a  rms  of  2.67  cm.  However,  the  agreement  in  the  analysis 
reported  here,  of  the  mean,  the  variance,  and  the  geographical  dependence 
of  the  SMMR  and  FNOC  data  would  support  confidence  in  the  SMMR.  This  issue 
can  be  resolved  by  investigation  of  additional  available  radiosonde  data. 

Second,  the  analysis  is  limited  to  the  Western  North  Atlantic  where 
the  FNOC  data  would  be  the  best  available.  In  most  other  oceanic  regions 
one  would  expect  poorer  quality  of  the  FNOC  hindcasts.  However,  since  the 
flight  of  Seasat-1  in  1978  the  general  quality  of  the  FNOC  hindcast  can  be 
expected  to  improve.  This  can  be  discussed  further  by  analyzing  the  re¬ 
maining  Seasat-1  data  and  a  general  examination  of  the  FNOC  data  products. 

In  particular  the  geographical  mean  wet  troposphere  correction  should  be 
obtained  for  each  of  the  world's  oceans. 

Third,  and  finally,  the  data  was  taken  in  the  three  month  period 
(July  -  October  1978).  During  this  summer  time  interval  (the  dog  days  or 
summer  doldrums),  one  would  expect  smoother  changes  than  one  would  during 
the  winter.  Therefore  the  variability  present  in  this  data  might  be  more 


benign  than  usual.  Furthermore,  the  short  time  interval  precludes  investi¬ 
gation  of  seasonal  variations  in  the  mean  tropospheric  correction  and  its 
variance.  This  seasonal  variation  must  exist  and  be  important.  Seasonal 
variation  of  the  mean  geographical  variation  can  be  investigated  using 
FNOC  data. 

Within  these  caveats,  one  can  conclude  that  using  the  FNOC  derived  wet 
troposphere  correction  will,  in  a  statistical  sense,  contribute  an  error  of 
(3.94  cm)**2  for  a  Seasat-1  type  of  data  span.  With  a  full  two  years, 
which  could  provide  14  points  in  each  15'  x  1 5 1  equatorial  square  after  two 
years,  one  can  anticipate  a  troposphere  error  of  (2.25  cm)**2. 

To  verify  and  insure  this  encouraging  result,  a  number  of  steps  should 
be  taken: 

1)  The  full  set  of  Seasat  SMMR  and  FNOC  data  should  be  analyzed  to 
obtain  an  initial  estimate  of  the  regional  distribution  of  the  mean  wet 
troposphere  correction. 

2)  The  accuracy  of  the  FNOC  data  in  other  oceans  should  be  examined 
and  regions  of  deficient  coverage  studied. 

3)  Regional  and  seasonal  models  of  the  wet  tropospheric  correction 
should  be  developed  for  parts  of  the  ocean  where  the  FNOC  data  does  not 
give  satisfactory  agreement. 

4)  The  Seasat  SMMR  calibration  can  be  redetermined  following  the  pro¬ 
cedure  of  Tapley  et  al  (1982).  However,  it  is  not  anticipated  that  signi¬ 
ficant  changes  in  the  calibration  will  occur,  though  a  more  reliable  (in 
the  sense  of  more  samples)  estimate  of  the  accuracy  of  the  SMMR  calibra¬ 
tion  will  be  achieved. 


SUMMARY 


We  have  investigated  the  accuracy  which  can  be  expected  using  Geosat-A 
satellite  to  sea  surface  radar  altimeter  data  for  determination  of  the  mean 
sea  surface  in  terms  of  the  ultimate  objective:  viz,  determining  the  marine 
geoid  and  gravity  at  sea.  This  view  allows  one  to  concentrate  on  regional, 
intermediate  wavelength,  time  invariant  error  sources.  This  is  important 
because  some  error  sources  that  have  long  wavelength  errors  along  a  track 
(e.g..  Orbit  errors)  have  random  (white)  spectra  when  averaged  over  space, 
i.e.,  when  viewed  across  track-to-track.  Other  error  sources  have  long 
wavelength  but  move  geographically,  and  have  a  similar  random  contribution 
when  averaged  over  time,  i.e.,  when  viewed  along  col  inear  tracks.  The  bene¬ 
fits  of  averaging  do  not  eliminate  the  desirability  of  modeling  all  effects 
with  utmost  care.  Clearly  one  would  not  want  to  rely  on  a  statistical  im¬ 
provement  if  not  necessary.  Listed  in  Table  9  are  the  Altimeter  Data  error 
sources  discussed.  We  have  already  discussed  the  errors  in  computing 
gravity  anomalies  as  a  function  of  geoid  height  (Sections  II,  IV,  V,  VI; 
and  Gaposchkin,  1983). 

We  assume  that  clean  data  is  available,  i.e.,  that  instrument  glitches, 
saturation,  specular  returns,  etc.,  are  already  removed.  Furthermore,  we 
assume  that  sensor  bias  depending  on  significant  wave  height,  and  attitude 
are  properly  calibrated.  In  this  case,  the  two  large  error  sources  in 
Table  6  are  the  orbit  error,  and  the  mesoscale  eddy  field.  Note  that  we 
believe  that  a  third  error  source  of  concern,  that  due  to  the  wet  tropo¬ 
sphere  correction,  does  not  seem  to  be  a  problem. 


The  orbit  error  has  been  discussed  above.  This  presents  a  major  chal¬ 
lenge.  No  doubt  by  application  of  methods  described  above,  significant 
progress  can  be  expected.  Nevertheless,  it  would  seem  unwise  to  assume  an 
orbit  error  much  better  than  10cm  along  any  given  track.  If  this  is 
achieved,  then  the  track-to-track  averaging  could  improve  this  by  about  a 
factor  of  two  to  three.  Therefore  we  can  hope  for  a  random  component  to  the 
orbit  error  of  3  to  5  cm. 

The  ocean  tide  models  are  quite  good  in  regions  where  tide  gauge  data 
is  available.  The  10  cm  uncertainty  is  really  conservative.  This  is  all 
one  can  say  in  the  absence  of  additional  information.  However,  the  tidal 
signal  is  time  variable  at  a  given  location.  Therefore,  even  if  there  was 
no  progress  in  tide  modeling,  which  is  unlikely,  the  natural  time  averaging 
for  an  individual  region  over  the  lifetime  of  Geosat-A  would  reduce  this 
error  source  to  less  than  2cm.  Of  course  there  will  certainly  be  local 
regions  with  large  tidal  signals  that  will  not  be  modeled  this  well. 

The  main  obstacle  to  complete  use  of  satellite  altimeter  data  to  deter¬ 
mination  of  the  marine  geoid  is  oceanographic  phenomena  that  cause  the  sea 
surface  to  differ  from  the  geoid.  Such  phenomena  have  amplitudes  of  a 
meter,  temporal  frequencies  of  days  to  decades,  and  spatial  frequencies  of 
50km  to  ocean  basin  scale,  10000km.  The  mean  circulation  can  be  eliminated 
by  avoiding  the  geographical  regions  of  the  strong  Western  Boundary  Currents 
such  as  the  Gulf  Stream  or  Kurshio.  The  ubiquitious  Mesoscale  Eddy  field 
can  be  approched  in  two  ways.  First,  by  using  the  natural  2  to  5  month 
time  constant,  averaging  over  the  full  two  year  Geosat-A  lifetime,  one  can 
reduce  this  100  cm  time  variable  signal  to  between  8  and  25  cm  rms.  Second, 
by  using  col  inear  tracks  and  simple  models  of  an  eddy  one  can  map  the  pro- 
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gress  of  the  eddy  field.  With  an  eddy  model,  matched  filters  could  be  suc¬ 
cessful  in  a  way  similar  to  the  study  of  seamounts  with  Seasat-1  data.  If 
successful,  this  model  can  be  subtracted  from  the  altimeter  data  before  it 
is  averaged  for  a  mean  sea  surface. 

Principal  Recommendations 

1.  Establish  the  global  geographical  and  seasonal  variation  of  the 
wet  troposphere  water  vapor  correction  using  the  Seasat-1  data  and  a  longer 
time  base  of  data  from  FNOC. 

2.  Develop  models  for  Mesoscale  Eddies  and  use  Seasat-1  as  a  test  for 
technique  development  to  map  the  evolution  of  the  eddy  field. 
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TABLE  9 


ERROR  BUDGET  FOR  GEOSAT-A 


Error  of  Single 
Measurement 

Phenomena  _ (cm) _ 


Altimeter  Noise 

2.5 

Sea  State  Bias  (2%  hl/3) 

2.  ■ 

Troposphere  dry 

0.7 

Troposphere  wet* 

5.0 

Ionosphere 

3.0 

Orbit  error** 

10.0 

Inverted  Barometer 

3.0 

Solid  Body  Tide 

2.0 

Ocean  Tide 

2. 

Mean  Ocean  Circulation*** 

100. 

Eddy  Field 

100. 

♦Assumes  FNOC  derived  correction. 

♦♦Assumes  Crossing  Arc  Orbit  Rectification. 
♦♦♦Only  Near  Western  Boundary  Currents. 


Error  Averaged 
Over  Mission 
(cm) 


2.2 

3.  -  5. 


2.0 

100. 

8.  -  25. 
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